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Abstract 

We derive source integrals for multipole moments that describe the 
behaviour of static and axially symmetric spacetimes close to spatial 
infinity. We assume isolated non-singular sources but will not restrict 
the matter content otherwise. Some future applications of these source 
integrals of the asymptotic multipole moments are outlined as well. 



1 Introduction 

In experiments that measure general relativistic effects, some parameters 
characterizing the spacetime have to be determined. The multipole mo- 
ments are one set of such parameters. They are measured in the exterior 
region of astrophysical objects like neutron stars or galaxies but also planets 
and describe the gravitational field near spatial infinity. They will be called 
here asymptotic multipole moments (AMM). The bending of light and the 
gravitational lensing proofed particularly useful for their measurement, see, 
e.g., [131 [HI [6] 8.nd references therein. But what information can be gathered 
about the matter distribution and the metric in its interior by their mea- 
surement? What does it mean to measure a quadrupole moment of a certain 
amount? In Newtonian theory, the answers are provided by the source in- 
tegrals of the AMM. These integrals determine the asymptotic multipole 
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moments by an integration over the mass density. In general relativity, sim- 
ilar expressions for the AMM are only known in approximations to general 
relativity, see e.g. [12] , Here we will present source integrals of the AMM 
for static and axially symmetric spacetimes in full general relativity. At the 
same time these provide quasi-local definitions of all asymptotic multipole 
moments. 

2 Preliminaries 

In this section, we shortly review several concepts necessary in the derivation 
of the source integrals. We use throughout this article geometric units G = 
c = 1 and the signature of the metric is (—,+,+,+). 

2.1 The line element and the field equations 

We concentrate on axially symmetric and static spacetimes of the Weyl form, 
i.e., 

ds' = e^^-^^ {dp" + d^) + H^2^-2^ V _ ^2U^^2^ 

We do not restrict the type of matter except in that the line element ([1]) can 
be introduced, see [TT]. The metric functions and W can be expressed 
by the timelike Killing vector = and the Killing vector of the axial 
symmetry ry" = (Si^)" 

Let us choose a sphere <So of finite radius r = Rq {p = rs'mO, ( = 
rcos9) that covers the entire matter distribution, cf. Fig. [TJ Outside of 5o, 
canonical Weyl coordinates {W = p) are introduced by virtue of one of the 
vacuum field equations. This allows still a shift in the C~coordinate, which 
enables us later to move the origin with respect to which the AMM are 
measured. The vacuum field equations read in canonical Weyl coordinates 

AU = 0, k^ = 2pUpU^, kp = p{{Upf-{U^f), (3) 

+ - ^ + ] . The function k is determined via a line 
integration, cf. the last two equations of ([3]), once U is known. Hence, only 
a Laplace equation for U remains to be solved in practise. 
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Figure 1: An example of the physical situations discussed here: The surfaces 

(2) (2) 

of the individual matter components are denoted by SI = dV^ with i>l 
and their respective volumes by V^-^^. The surface 5q^^ describes a circle 
with sufficiently large but finite radius enclosing all matter components, cf. 
Sec. EH 

2.2 The physical setting 

We depict in Fig. [1] an example of a physical situation that will be covered 
by the subsequent considerations. The relevant surfaces and volumes are 
defined there as well. For simplicity, we allow only non-singular sources. 
However, such can be incorporated into the formalism as we showed in [5]. 
The 3-dimensional projection of the matter region into an hypersurface of 
constant Killing-time t is obtained by an rotation around the axis in Fig. 
[TJ In this way, the quantities S^"^ and V^^^^ are defined starting from 
^^"^ H'^'*' respectively. 

2.3 The multipole moments 

For asymptotically fiat and static spacetimes a geometric definition of AMM 
was given by Geroch in [3]. This definition was generalized and applied by 
many authors, see the reviews [121 ED] aiid references therein. In the axially 
symmetric case with the line element ([T]), Geroch's multipole moments My. 
can be obtained by an expansion of U along the symmetry axis in |C|~^ '■ 
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The follow uniquely from Weyl's multipole moments C/'-^^ and vice versa 
as was shown in [2]. Therefore, we consider only the t/^^-* here. 



2.4 The inverse scattering technique 

Lastly, we shortly review the inverse scattering technique (1ST), see e.g. f7] 
for a recent account. Even though the Laplace equation is linear and the use 
of the 1ST seems artificial, the 1ST proves nonetheless beneficial, because 
it is easily generalizable to the non-linear case of the Ernst equation. This 
equation is of special interest in relativistic astrophysics, since it describes 
the exterior of rotating stars. The starting point of the 1ST in the present 
setting, i.e. the linear problem of the Laplace equation, is given by 

a , = (1 + A)C/,,(T, a , = (l + ^) U-^a, (5) 



where z = p + i^, the spectral parameter A = y G C and a bar 

denotes complex conjugation. The complex valued function a depends on 
z, z and A. The integr ability condition of Eqs. ^ is the Laplace equation 
for U . Therefore, having a solution a of Eq. ([5]) yields also a solution U 
of the Laplace equation and vice versa. The main technical steps of the 
1ST as described in are to integrate ([5]) along Aj^ , along a circle with 
sufficiently large radius and along a compact curve connecting A\ with 

A_ . This scheme can be carried out partially and we quote only the results 
(simplified to the static case), which are relevant for us, from |7j: 

(0, C) G ^+ : a{\ = +l,p = 0, C) = F(if)e2^(''=0'^), 

a(A = -l,/5 = 0,0 = 1, 
(0,0 G : a (A = +1,/; = 0,0 = e2^(''=0'«, 

a{\ = -l,p = 0,O = F{K). 



(6a) 



The function F : C ^ C is given for G M with {p = OX = K) eA^ by 

(2) 

The integration along Sq does not enter the derivation of these formulas 
and it forms the crucial part of our considerations in the next section. 
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3 Source integrals of Weyl's multipole moments 



The derivation of the source integrals consists of several steps. First, the 

(2) 

AMM are expressed as line integrals along 5q . This is the most important 
step, because it makes it possible to determine the AMM quasi-locally. Then 
these integrals will be rewritten in a coordinate independent form as surface 

(3) 

integrals over Sq by virtue of the axial symmetry. Subsequently, Stokes' 
theorem is used to rewrite these as volume integrals over V(°). In the final 
step, it is shown that the contributions in the vacuum regions vanish. Thus, 
the steps from before can be retraced to obtain the contributions in source 
integral form of each individual matter component. We suppress the details 
of these derivations and show only the crucial steps. 

(2) 

The linear problem ([5]) is well-defined along Sk and reads: 



(7) 



where U^s and U^n are the tangential and the (outward pointing) normal 

(2) (2) 

derivative of U along Sq with respect to a parametrisation [sn, ss] — ?• 5q . 
The indices N and S refer to the values of a parameter or function at the 

(2) 

"north" and "south" pole of 5q , i.e., to the intersection points (p = 0, C = 

(2) 

Cn/s) of iSg and the symmetry axis. Eq. ([7]) is easily integrated using the 
boundary values from Eq. ([6]): 

^(0,i^)= ^^~^^ +l J ((X'^ + x)Us + i{\-'-X)Un)ds. (8) 

If we expand this equation in |-ftr|~^, we obtain expression for Weyl's mul- 
tipole moments in terms of a line integration. Let us introduce the ab- 
breviations N^^ and N^^ for the expansion coefficients (A~"^ + A)*^^^ and 
i(A~"'^ — A)'-'"^ to order r + 1, respectively. After a lengthy but straightfor- 
ward calculation they evaluate to 

4fc(A;!)2(r-2A;)! ' 

(9) 



zL> 4'=(A;!)2(r-2A:- l)!(2A: + 2) 
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The r = — 1 order in l-R'l^^ of Eq. ([8|) is satisfied trivially and will not 
be considered subsequently. The orders r > of Eq. ([8]) yield the desired 
quasi-local definitions of Weyl's multipole moments: 

U^''^ = \j + ivlrV,,) d5f , (10) 

where U^s and U^fi are the tangential and normal derivatives along Sq with 
respect to the unit tangent vector and the unit normal vector, which are 

(2) (2) 

defined with the induced metric on 5q . d^Q denotes the proper distance 
along 5o^^ . The functions N^^^ and U are to be read as functions along s'^^ , 
i.e. as functions of {p{s) ^ C,{s)) . 

To make the coordinate independence apparent, we express p and C by 
scalars build from the Killing vectors. Firstly, observe that Eq. ([2|) holds 
everywhere and that in vacuum we have W = p. Additionally, the 1-form 

Za = eat3y5W'^W-^rjH^ (H) 

is well-defined and hypersurface orthogonal everywhere as well as exact in 
the vacuum region. Hence, there exist a potential Z and an integrating factor 
X such that Z^a = XZa, where X = 1 in the exterior of In the vacuum 
region and in canonical Weyl coordinates, we find Z = + const. Since we 
can shift the ^-coordinate freely, we can drop the constant of integration, 
which specifies the origin with respect to which the AMM are measured. 
Thus, W and Z coincide with p and C in the vacuum region and can be 
used as their continuation into the interior of the matter. This choice is not 
unique and other continuations are possible, although they do not alter the 
values of the source integrals, which we present below. 

(2) 

Using W and Z along Sq instead of p and respectively, we can rewrite 
Eqs. (flOl) as surface integrals: 

U^'^ = il"^ {n^-^U, - N^:^^Z,,U + d5f , (12) 

A3) 
^0 

An integration by parts, the axially symmetry and the vacuum field equa- 
tions are necessary for this step. 
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Using Stokes' theorem and the field equations we obtain 



N^""^ (r) (W 
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;(3) 



W 



dV, 



(3) 



w 



dV) 



(3) 



(13) 



The dvf ^ are the proper volume element of V^^^^ and a semicolon denotes the 
covariant derivative with respect to the line element ([1]). The last equality 
is due to Einstein's equations, which imply that the integrand vanishes in 
vacuum. The integrals (|13p are the desired source integrals. They determine 
the AMM from the geometry inside the matter regions alone. Of course, 
Stokes' theorem can again be used to rewrite the source integrals as surface 

(3) 

integrals over 5^ of the respective matter component. In turn, these can 
be reformulated as line integrals, cf. Sec. HI That the contributions of the 

(3) 

individual matter components to the asymptotic multipole moments 
superpose linearly is due to the choice of Weyl's multipole moments. If we 
employ the method from [2] to calculate Geroch's multipole moments M^. 
from Weyl's multipole moments U^'''\ we obtain a mixing of the contributions 

(k) 

U- of the individual matter components with k < r in the M^. This 
is already apparent for the quadrupole moment M2, which depends non- 
linearly on U^^'^: 



M2 = - 

3 



(14) 



The Geroch mass Mo equals [/(°) and is given by the (negative) Komar 
integral. This follows also from Eq. (jl2p with r = 0. 
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4 Applications 

We conclude the paper discussing one possible application of the source 
integrals (jl2p . Assume a matter distribution is given, where the metric is 
known in the interior or the Dirichlet and the Neumann data for U are known 
at the surface. Even then it is far from trivial (at least in the stationary 
case, see [IJ) to obtain a global asymptotically flat solution, if it exists. 
The source integrals for the AMM provide a tool to solve this task. As 
a simple example serves here the case of static dust without any surface 
distributions. In Weyl coordinates (not necessarily canonical) the energy 
momentum tensor is given by 



The contracted Bianchi identities imply U = const, in the interior and, thus, 

(3) 

the gradient of U vanishes at in all coordinates. Using the line integrals 
for Weyl's AMM, which follow from Eq. (fT3l) . we get: 



Thus, all AMM vanish and the spacetime is flat in the exterior. This con- 
tradicts the presence of a dust distribution with positive mass density. Of 
course, this result is already known and more general non-existence results 
for dust including the rotating case can be found in [H |9] and references 
therein. Although the non-existence is proved here, this example shows in a 
concise way how the source integrals can be applied in more difficult phys- 
ical situations like rotating stars. This and other applications, e.g. to tidal 
distortions of black holes, will be investigated in future work. 
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